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ABSTRACT 

A putative temporally varying circulation-free magnetic-field configuration is inferred 
in an equatorial segment of the solar convection zone from the helioseismologically 
inferred angular- velocity variation, assuming that the predominant dynamics is angu- 
lar acceleration produced by the azimuthal Maxwell stress exerted by a field whose 
surface values are consistent with photospheric line-of-sight measurements. 

Key words: Sun: dynamics; Sun: helioseismology; Sun; interior; Sun: magnetic field; 
Sun: rotation 
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1 INTRODUCTION 

This paper reports a first step towards investigating directly 
from helioseismological inference of angular-velocity varia- 
tions the magnetohydrodynamical force balance in the solar 
convection zone. The balance is predominantly azimuthal, 
and we presume that rotational acceleration is caused by 
Maxwell stresses. We ignore all other forces. 

Our long-term goal is to infer from seismic and surface- 
magnetic-field data the distribution of the axially symmetric 
component of the magnetic field throughout the convection 
zone, and possibly beneath, together with the meridional 
circulation associated with it. Our intention is to accom- 
plish that task by an inverse dynamical argument, starting 
from what one might regard as an outcome of the dynamics, 
namely the temporally varying angular velocity, which has 
already been determined seismologically. The task is thus to 
find a magnetic-field configuration and a circulation pattern 
that both satisfy the magnetic induction equation and pro- 
duce a Maxwell stress which, together with the transport of 
angular momentum by both the meridional fiow and puta- 
tive Reynolds stresses, induce the observed angular velocity 
as a consistent consequence of the azimuthal component of 
the (angular) momentum equation. 

Buoyancy-driven motion in the convection zone redis- 
tributes angular momentum and magnetic field, the larger- 
scale motion - what is normally termed rotation, zonal 
fiow and meridional circulation - via laminar advection, the 
smaller scales via anisotropic Reynolds stresses and turbu- 
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lent diffusion. Usual (forward) theoretical investigations are 
aimed at reproducing that motion and its consequences, such 
as the temporal variation in the angular velocity, which is 
then to be compared with observation. To that end it is 
necessary to solve the full coupled momentum and magnetic 
induction equations. Our intention here is different. We ar- 
gue inversely: to infer from the observed variation in the 
angular velocity what force is required to drive it. In this 
first foray we assume axisymmetry; in that case only the az- 
imuthal component of the momentum equation is required 
for inferring the pertinent components of the stress. It is 
not necessary to consider the meridional components of the 
momentum equation; in principle they can be used subse- 
quently to infer from our results the forces that control the 
fiow, and thereby provide a valuable diagnostic for forward 
dynamical calculations. That lies outside the domain of our 
present endeavour, so we do not discuss it here. Nor do we 
discuss the underlying mechanisms that maintain the varia- 
tion observed. 

The problem thereby presented to us is highly nonlin- 
ear. Our hope is that it can be solved by iteration. That 
will not be straightforward, partly because it is necessary 
first to obtain a trial that is close enough to the solution to 
render the procedure convergent. We need also to discover 
an appropriate iterative procedure. To this end we first at- 
tempt to find a solution with no meridional fiow. And we 
go no farther than that in this preliminary exercise. That 
task is also nonlinear, although apparently less so. We have 
succeeded so far in solving it by iterative back-substitution 
in an equatorial region occupying latitudes up to 45°. 

The primary data for our investigation are the tempo- 
ral variations over the past one-and-a-half solar cycles of the 
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angular velocity Q in the solar interior, obtained from rota- 
tional splitting of seismic modes produced by advection and 
Coriolis precession. The most prominent feature of the varia- 
tion is a series of latitudinal bands of faster- and slower-than- 
avcrage rotation, which match the so-called torsional oscilla- 
tions observed by Howard and LaBonte (1980) in the photo- 
sphere. At low latitudes, the bands appear to drift towards 
the equator; at high latitudes they migrate polewards (Howe 
et al. 2000; Antia & Basu 2001). They persist throughout 
the convection zone, and possibly extend into the radiative 
region immediately beneath (Vorontsov et al. 2002; Basu & 
Antia 2003; Howe et al. 2005; Antia, Basu and Chitre 2008). 
The pattern seems to be well correlated with the magnetic 
butterfly diagram (Antia, Basu and Chitre 2008; Sivaraman 
et al. 2008), and it is commonly believed to result from the 
interaction between differential rotation and magnetic field 
(e.g., Covas et al. 2000; Rempel 2006). Thus the properties 
of zonal flows provide a strong constraint on the dynamics 
of the interior of the Sun. 

The frequencies w of the seismic oscillations are also af- 
fected directly by the magnetic field: modes of like order n 
and degree I are split with respect to azimuthal order m. 
The frequency splitting is very small, so it would be diffi- 
cult to infer the structure of the field even if it were known 
that it is a magnetic field that is the cause of the splitting. 
Like the effect of any other purely structural asphericity, the 
splitting that is induced is an even function of m, because 
the Lorentz force cannot distinguish between cast and west. 
Indeed, Zweibel and Gough (1995) have argued that one can 
never distinguish a magnetic field from other structural as- 
phericity from seismic frequencies alone. Rotation, on the 
other hand, is axial, and docs so distinguisfi; it gives rise to 
frequency perturbations with an odd component - the term 
first-order in Q/u) is purely odd in m, and is used to infer 
Q{r,9,i); once that has been determined, the second-order 
perturbations can be calculated and subtracted from the ob- 
served even splitting to yield a residual contribution which 
is due to the remaining component of structural asphericity, 
and which could be the direct result of a magnetic field. If 
one were then to assume that the resultant even contribu- 
tion is due solely to magnetism it may be possible to make 
some direct seismological inference about the structure of 
the field. 

The contribution to firequency splitting due to a large- 
scale magnetic field can be calculated using the formulation 
given by Cough & Thompson (1990). From that analysis, 
Antia, Chitre & Thompson (2000, 2003) have inferred the 
possibility of a magnetic field of about 70 kG in the vicinity 
of r = O.9-R0 (cf. Kosovichev, 1996), and an upper limit of 
300 kC on a possible field near the base of the convection 
zone. Chou & Sercbryanskiy (2005) used frequencies from 
the Michelson Doppler Imager (Scherrer et al. , 1995) to de- 
tect some temporal variation correlated with solar activity 
which could be duo to magnetic-field variation near the base 
of the convection zone of intensity several tens kC, and Bald- 
ner et al. (2009) have tried to infer magnetic-field variation 
in the near-surface regions from seismic frequencies alone. 

Here we address the rotational dynamics directly. We 
enquire what force is required to induce the observed vari- 
ation in the sun's angular velocity, under the assumption 
that the driving force is solely magnetic, an assumption that 
was made in early studies of the dynamics of the torsional 



oscillations (Schiissler, 1981; Yoshimura, 1981), and we de- 
duce the magnetic field that produces it. In this tentative 
foray we ignore viscous and Ohmic diffusion, be they lam- 
inar or turbulent. We also ignore meridional advection and 
its contribution to the distortion of the magnetic field. We 
recognize that this highly simplifying assumption does not 
take into account all the pertinent dynamical processes that 
might account for the variation, but we present it as the first 
step towards obtaining a dynamical description of the solar 
cycle that is consistent with observation. 

An objective superficially similar to ours appears to 
have been the motivation of a recent study by Lanza (2007). 
Lanza does not adopt the hclioscismologically inferred an- 
gular velocity, but instead assumes that the angular velocity 
in the convection zone results simply by pure diffusion from 
the angular velocity observed at the solar surface. Therefore 
his analysis bears little resemblance to ours. In his computa- 
tions ho ignores meridional flow, as do we. Moreover, Lanza 
assumes that the principal component of the Maxwell stress 
that drives the flow is BrB^, and he ignores BeB^; we take 
both into account, and we find that in most places it is the 
latter that dominates. 

Our procedure is to force a balance between the az- 
imuthal component of the Maxwell stress and the observed 
angular acceleration. To that end, we use the seismologi- 
cally inferred angular velocity i}{r,6, t) that was used in a 
recent study of the temporal variations in the Sun's angu- 
lar momentum and rotational kinetic energy (Antia, Chitre 
& Gough 2008). In relating the Maxwell stress to the mag- 
netic field, B, we confine attention to low latitudes where 
we presume (initially) that the latitudinal component of the 
field dominates over the radial component. That permits an 
iterative procedure starting from a state in which the radial 
field is ignored, ending with a state in which all field com- 
ponents are consistently accounted for. At high latitudes it 
is the radial component that dominates; in principle a cor- 
responding procedure might also be applied, but there the 
angular velocity is less well determined and therefore so is 
the Maxwell stress; moreover, the mathematical structure of 
the iterative procedure is different. We suspect that conver- 
gence is less secure. 

Since angular velocity is the only property of the 
zonal flow that we have observed, we must characterize 
the Maxwell stress with an effective axisymmetric magnetic 
field. Moreover, because only the north-south symmetric 
component of O is available from the inversions of seismic 
frequency data, we adopt a corresponding symmetry in B, 
which, in accord with observations of the large-scale dipo- 
lar component in and above the photoshere, we presume to 
change polarity between sunspot cycles. In this investigation 
we do not attempt to satisfy the latitudinal component of 
the momentum equation. Instead, we take the observed an- 
gular velocity as our starting point. The component of the 
Maxwell stress that drives that flow also drives a meridional 
circulation, which advects both magnetic field and angular 
momentum (moderated by diffusion), thereby complicating 
the situation considerably. We have not succeeded in de- 
termining a diffusion-free flow that satisfies the consequent 
dynamical balance; we recognize that one may not even ex- 
ist. One can easily calculate a flow that arises from the field 
we have calculated by ignoring advection and diffusion, but 
subsequent iteration to accommodate advection based on 
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straightforward back-substitution, perhaps unsurprisingly, 
does not converge. To make further progress therefore de- 
mands a more sophisticated approach. 

We relate the Maxwell stress to the latitudinal and az- 
imuthal components of the magnetic field B via the cor- 
responding components of the magnetic induction equation 
from which, subject to our assumptions, the horizontal com- 
ponent of B is determined in terms of an integral of SI. 
The radial component of B is then determined from the 
divergence equation V ■ B = 0. Completing the solution re- 
quires the imposition of boundary conditions; we use the az- 
imuthally averaged observed line-of-sight component of the 
magnetic field in the solar atmosphere. 



2 FORMULATION 

We ignore the meridional flow velocity together with pertur- 
bations to the hydrostatic stratification induced by that flow 
and by the angular velocity, and indeed also the magnetic 
field itself. Thus, with respect to spherical polar co-ordinates 
(r, 6, (j>) referred to the axis of rotation, the azimuthal com- 
ponent of the diffusion-free momentum equation for an ax- 
isymmetric configuration is 



dv^ _ 1 d 
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where t is time, v^f, = rfl sinO is the velocity of rotation, and 

r = BeB^., A = (2) 

in which B — {Br,Bg,B^). The corresponding component 
of the diffusion-free magnetic induction equation is 
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We confine attention to latitudes equatorward of 45° , where 
we expect IB^I <g \Be\, implying |A| ^ jr|. We assume also 
that the meridional field has the same symmetry as a dipole: 
Br{r, 71-9,1) = -Br{r,9,t) and Be{r,-K-9,t) = Be{r,9,t), 
implying that the azimuthal component, which is generated 
by stretching due to rotational shear with S7(r, tv — 9,t) — 
Q{r,9,t), satisfies B^(r,9,t) = —B^{r,n — 9,t), whence 
B^{r,T7 /2,t) = 0. Under these symmetries, equation ([T} can 
be integrated with respect to 9 to yield 

f.7r/2 
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which is then combined with equation Q and integrated to 
obtain B^. In carrying out that procedure we notice that 
Q{r, 9, t) may be written as f2o(r-, ^) + Sli (r, 9, t) with \Q,i \ <^ 
f2o|. Accordingly, it is expedient to use this separation for 
dQ,/d6 to obtain 
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in which to is some fiducial origin of time; H{r, 9) is a func- 
tion of integration, depending, of course, on the choice of 
to- For to we use the instant at which seismic data are 
available. In equations (6) and (8) we have neglected the 
temporally varying terms in dQ/d9 and dfl/dr to empha- 
size from where the major time dependence in the solution 
arises; the approximation modifies coefficients multiplying 
relatively small terms by no more than 0.3%, which makes 
an imperceptible difference to the solution at the resolution 
of Figures 1-3. Given B^, the other horizontal component 
of the field can be obtained as F/B^ from equation ([SJ, and 
then Br can be obtained by integrating V ■ B = 0. We point 
out that one might have expected an additional function of 
integration, namely M{r,t), added to the right-hand side of 
equation ([5} ; it is straightforward to show that this function 
leads to a component of the magnetic field that does not 
satisfy the required symmetries, and must therefore be null. 

We anticipate that at any given location B^ changes 
sign at some instant during the solar cycle. That can occur 
only when _B| = 0. Therefore, we choose the function H of 
integration to render it possible: 



(9) 



The sign of B^ is chosen to match the sign of the observed 
field near the surface; Bo changes sign at other instants at 
which r vanishes. 

We introduce a flux function, x, which satisfies 

(S.-S«)-^(ig--§^), (10) 
r sm9 \r o9 or / 

which can be integrated to 

X{r,9,t) = sm6 f r' Be dr' + x{R,9,t) , (11) 

J r 

in which the second term is determined from the longitudi- 
nally averaged line-of-sight magnetic field B^^{9,t) observed 
at the surface r — R of the Sun: 



r/2 



x{R,9,t) = R^ / {Be cos 9' - B\^) d9' . 



(12) 



These equations have been solved by iterative back- 
substitution: expressions ©-0 and ® for F and were 
evaluated first with A = and A = 0, from which a trial Be, 
and subsequently x and Br, were obtained from equations 
(HI, (fTTj) and (fT^ . and ^TU^. New values of A and A were 
then evaluated, and the whole procedure repeated till the 
solution converged. 



3 SEISMIC DATA AND SURFACE MAGNETIC 
FIELD 

We have used GONG frequency-splitting data (Hill et 
al. 1996) collected over a period from 1995 to 2011 (GONG 
months 1-165). Each data set covers three overlapping 
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GONG months (of 36 consecutive days each), with mid- 
points separated by 1 GONG month. From each set the 
angular velocity Q, was inferred as a function of radius and 
latitude by 2-d RLS data fitting, implemented in the man- 
ner of Antia, Basu and Chitre (1998). The results were then 
smoothed at each radius and latitude by fitting them, by 
least squares, to the oscillatory form: 

2 

n(r, 9, t) = (n) + ^ ttfe sin(A;a;ot + (j)k) =■ + 50. , (13) 
fc=i 

in which the angular brackets denote average over time, and 
Wo is the frequency corresponding to a period of 11. 7y, the 
inferred length of solar cycle 23 (Antia & Basu 2010); thus 
(fj) represents Oo and 50, represents f2i in equations (6) 
and (8). The temporal expansion was taken to second order 
in order to accommodate quadratic nonlinearities. Note that 
(Q), Ofc and 4>k are functions of r and 0. The smoothed angu- 
lar velocity and its temporal derivatives were used to deduce 
the temporally varying magnetic field using the procedure 
described in §2. 

Full-disc magnetograms from the Wilcox Solar Obser- 
vatory have long provided valuable data for global magnetic 
studies (e.g. Duvall et al., 1979). Here we use data obtained 
by Lo, Hoeksema and Scherrer (2010), kindly provided by 
T. Hoeksema in the form of a longitudinally averaged line- 
of-sight magnetic field B\\{0,t) — the zeroth moment of the 
magnetograms with respect to longitude - which constitutes 
our surface boundary condition. 



4 RESULTS 

The starting point of our analysis is the observed temporally 
varying component 5Q.{r, 9,t) of the angular velocity. It is 
depicted in Figure 1 in the equatorial sector of a meridional 
quadrant within which our analysis has been carried out, at 
epochs separated, in the main, by two years. At amy instant 
there are typically three principal regions, with respect to 
latitude, of coherent variation; they are of alternating sense, 
and they extend over the entire vertical range of the con- 
vection zone, although the boundaries are not vertical: the 
pattern gradually distorts. However, the predominant vari- 
ation is a migration towards the equator at a rate which 
causes at (almost) every fixed location an oscillation with 
the period of the sunspot cycle, as has been observed at 
the the surface of the Sun (Howard and LaBonte, 1980) and 
in the subphotospheric layers of the convection zone (e.g. 
Vorontsov et al. 2002). The sign of 50 changes rather later 
than the sunspot maximum and sunspot minima. The max- 
imum amplitude of 50 is a little above 1.2 uHz. Note that 
the Lorentz force density, L^, also shows temporal varia- 
tion, and has typically two main regions in latitude extend- 
ing over the entire vertical range of the convection zone. Its 
maximum magnitude occurs around O.SRq. 

The (iominant Maxwell stress that drives the flow is 
proportional to F = BgB^ (see equation (1)). It is displayed 
in Figure 2, alongside the other pertinent component, A. 
There is little obvious large-scale migration in the convec- 
tion zone, unlike the variation of 50 from which F was de- 
rived; rather, simply an alternation in sign, occurring a lit- 
tle before the sign change in 50. However, the sign change 



does actually begin as an equatorward migration of a tongue 
of oppositely signed stress centred at r/R ~ 0.8, emerging 
into our region of consideration from the highest latitude in 
1996, 2002 and again in 2008 (a panel displaying the last of 
these is included in Figure 2; it is evident in Figure 3 that 
it stems principally from a similar configuration ol B^), and 
which rapidly expands to occupy much of the convection 
zone. During this process, the stress in the convection zone 
typically suffers two sign changes with respect to radius, 
particularly at the higher latitudes. Note that the contours 
of A = BrB^ are substantially more widely separated than 
those oiV = BeB^, indicating that on the whole the radial 
component of the meridional field has a lesser influence than 
the horizontal component, as we anticipated. The stress dis- 
tribution in the tachocline beneath the convection zone is 
more complicated, possibly because advection of the mag- 
netic field by a meridional circulation - not incorporated 
here - is important there; broadly speaking, the stress is in 
the opposite sense to that in the convection zone immedi- 
ately above, as is to be expected. 

Contours of constant flux function, x-i which map the 
poloidal magnetic field, and contours on constant B^, are 
depicted in Figure 3. The poloidal field is predominantly 
dipolar near sunspot maximum, around 2000, and predom- 
inantly quadrupolar, with higher multipoles also evident, 
near sunspot minimum, both near 1996 and throughout the 
latest broad minimum from about 2005 to the present. The 
dipolar component reverses sign during the declining phase 
of the cycle, somewhat earlier than sunspot minimum; it 
appears to have a characteristic period twice that of the 
sunspot cycle, approximately 22 y, as is observed to be 
the case at the surface. The azirrmthal field reverses sign 
throughout most of the convection zone near sunspot min- 
imum, and it too has a dominant period twice that of the 
sunspot cycle. Its maximum intensity, near the 2001 sunspot 
maximum, is about 4.5 kG; today it has attained the same 
value, and is still rising as the next maximum is approached. 
The field in the tachocline is more complicated: only near 
sunspot minimum is there a significant radial field, which 
then appears to be wound up by the shear to generate az- 
imuthal field; near sunspot maximum the tachocline field is 
predominantly horizontal. 

The temporal variation of the field is resolved better 
in Figure 4. There Be and B^ are plotted against latitude 
and time at various radii. Perhaps the most striking fea- 
ture is the upward migration of the pattern of Bg, which is 
most easily seen in the location and shapes of the null lines: 
the time lag behind the variation at the base of the con- 
vection zone is greatest near the equator, reaching a value 
of about 14y at r/R = 0.98; the lag is least at latitudes 
±30°, the latitudes at which sunspots first appear at the 
surface. In addition there is a temporal overtone in Bg near 
the equator, which appears to be related to the octupolar 
component of the poloidal field (cf. Figure 2). There is also 
an evident temporal overtone at mid latitudes near the base 
of the convection zone. The azimuthal component B^ of the 
field also migrates upwards, although only above r/R ~ 0.8; 
the greater migration speed near latitudes ±30° is perhaps 
somewhat more obvious than it is for Bg. The azimuthal 
component of the field changes sign near the surface around 
2006, near the beginning of the prolonged period of low ac- 
tivity from 2005 until about 2011. The reversal of Be near 
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Figure 1. Contours of the angular- velocity variation 5Q and the 
Lorentz force density (defined by equation (1) ), each averaged 
over three GONG months in the middles of the years indicated in 
the panels. Red contours are positive (i.e. they correspond to az- 
imuthal motion and force density in the direction of the mean 
rotational motion), blue negative. The continuous black circular 
segment marks the photosphere, the dashed segment the base of 
the convection zone. The contour spacings are 0.1 nHz in SQ/2tt, 
and 5 X 10~* gcm~^ s~^ in L^, except that the null contours are 
not dravirn. 



the surface also starts near the beginning of the prolonged 
sunspot minimum, at latitudes near ±30°, but occurs later 
at low latitudes, extending throughout the period of low 
activity. The greatest intensity, Be, max, of Be, was about 
lOOG, and occurred at the beginning of the period of inves- 
tigation near the the base of the convection zone; judging 
from the upward migration rate of the field, it is likely that 
the actual maximum had occurred about 2y before that. The 
radial component, Br, is considerably less intense, so Be, max 
estimates the maximum intensity of the total poloidal field. 
The intensity of the poloidal field appears to have attained 
a second maximum in 2011, with a value of only about 55G. 



5 DISCUSSION 

We have attempted to infer the magnetic-field configuration 
and its temporal variation in the solar convection zone by 
dynamical reasoning, starting from the seismically inferred 
angular velocity. Our analysis is based on the assumption 
that the Maxwell stress is wholly responsible for driving 
the time-varying axisymmetric component of the migrating 
bands of zonal fiow. We have ignored viscous and Ohmic 
diffusion; we have also ignored meridional fiow. 

It should be appreciated that the magnetic field in the 
Sun may be inhomogeneous on small scales. Therefore the 
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Figure 2. Contours of the components F = BgB^ and A = BrB^ 
of the Maxwell stress, each inferred over three GONG months in 
the middles of the years indicated. Red contours are positive, blue 
negative. The contour spacings are (100 G)'^ in both, except that 
the null contours are not drawn. 



field intensities that we have quoted are essentially rms val- 
ues characteristic of the spatially averaged Maxwell stress. 

Some features of our findings are in accord with com- 
mon expectation: an oscillatory, predominantly dipolar (at 
times of sunspot maximum), poloidal magnetic field with a 
22-y period, which reverses near sunspot minimum, reveal- 
ing a predominantly quadrupolar residual component near 
reversal; and a large-scale toroidal field, presumably gener- 
ated by the stretching of the poloidal field by the differential 
rotation, which also has a period of 22y. What is new is the 
seismological evidence that the reversal of the magnetic field 
migrates upwards from the base of the convection zone. 

The variation of the field that we find in the tachocline is 
more complicated. Comparison of the radial and azimuthal 
components of the field suggests that there is a signifi- 
cant meridional flow in the tachocline which oscillates with 
the cycle and which advects the fleld: equatorward during 
the rising phase of the cycle, poleward during the declin- 
ing phase. However, we must emphasize that a fully consis- 
tent calculation incorporating advection has not been car- 
ried out, so the hint must be taken with a pinch of salt. 

Our omission of meridional circulation is a serious de- 
ficiency in our calculations. We have assumed that the pre- 
dominant force causing the rotational acceleration arises 
from Maxwell stresses. But if a major component of the 
angular-momentum variation were to be due to advection, 
our results would prove to be inadequate. Is that actu- 
ally the situation? Our failure to obtain convergence by 
straightforward iterative back-substitution, starting from 
the circulation-free solution described in this paper, and ob- 
taining an estimate of the circulation at each iteration from 
the magnetic field configuration via the meridional compo- 
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Figure 3. Contours of constant flux function x, which indicate 
the direction of the poloidal component of the magnetic field, 
and the azimuthal component of the field, using the same 
conventions as in Figure 1: the direction of the poloidal field is 
clockwise on red contours of x, anticlockwise on blue. The contour 
spacing is 0.5 GRg in x ^i^d 500 G in B^, except that the null 
contours are not drawn. 
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Figure 4. In the left-hand panels are displayed contours of con- 
stant Bg inferred from three- month averages of SQ{r, 8, t) ob- 
tained by inverting all the (overlapping) GONG data sets and 
fitting to expression (13), plotted against epoch and latitude 
at the four radii indicated; in the right-hand panels are dis- 
played B^. The contour separations are IG, 2.5 G, 5G and 
5G in Bg, and 100 G, 250 G, 500 G and 500 G in at r = 
0.98, 0.90, 0.80 and 0.70 Rq, respectively; red represents posi- 
tive values, blue negative. The plots in the southern hemisphere 
were obtained from the assumed symmetry relations quoted in 
the text. As in the previous two figures, the null contours are not 
drawn. 



nents of the induction equation, suggests that it may be. 
But studies of simplified, yet dynamically consistent, dy- 
namo models by Rempel (2006) suggest that the influence of 
circulation on the overall structure of the field may not be se- 
vere. That supports our optimism that a relatively straight- 
forward dynamical inversion procedure might be within our 
grasp. What might be the necessary ingredients? 

First we intend to extend our circulation-free solution 
to the poles. Although we have found that it is possible to 
extend the procedure described in this paper to somewhat 
higher latitudes, it cannot be extended far, because the as- 
sumed low Br'.Be ratio adopted to justify the neglect of Br 
in computing our initial estimate of the field in the equa- 
torial regions is violated near the poles. Instead, we should 
design a procedure based on an initial state obtained by 
ignoring Bg; that leads to a problem with a different math- 
ematical structure. It is interesting to observe that in com- 
pleting the computation of the field structure in that case 
it is the first moment with respect to longitude of the the 
observed line-of-sight field that is needed for providing the 
magnetic boundary condition at the surface. That is math- 
ematically - although hardly physically - independent of 
the zeroth moment used in the equatorially concentrated 
analysis presented here. Therefore, the degree to which the 
polar and equatorial solutions match at mid-latitudes will 
provide some measure of the degree to which our underlying 
assumptions are valid. 

Next we shall attempt to incorporate the meridional 
flow V. That is to be accomplished by requiring that it ad- 
vects, against Ohmic diffusion, the magnetic field in such a 
manner as to yield the field B whose azimuthal stress, to- 
gether with angular-momentum advection by v, maintains 
the the observed angular-velocity variation 5il. It is deter- 
mined by the 9 component of the induction equation cou- 
pled with the continuity equation V.(pv) = 0. To complete 
the integration requires an additional boundary condition: 
namely, the value of the product BrVg at the surface, which 
can be obtained with the help of either the seismologically 
inferred meridional fiow near the surface (and extrapolated 
to the surface) or observations of small photospheric tracers. 

If eventually we succeed in solving the full problem, 
what observable predictions might there be for testing the 
outcome? We cannot use the meridional fiow observed in the 
photospheric layers because that will be needed as a bound- 
ary condition for completing the integration for the stream 
function. However, the fiow inferred well below the photo- 
sphere could be compared with direct helioseismological in- 
ferences, which are likely to have been extended to substan- 
tially greater depths than have been achieved to date. In par- 
ticular, it would be interesting to see whether we reproduce 
the double-cell structure suggested by recent time-distance 
analysis (J. Zhou, personal communication, 2012). It would 
be interesting also to see whether the inferred fiow matches 
the upwelling from the tachocline predicted by Spiegel and 
Zahn (1992) and Gough and Mclntyre (1998), although that 
component of the fiow might be too weak to have a notice- 
able impact on the fiow geometry above. In any case, our 
eventual results would provide a serious template against 
which to compare dynamical theories of the convection zone. 

Finally we recall that we have assumed axisymmetry 
on the large scale. To study asymmetric deviations would 
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require seismic inference of the noncixisymmetric zonal flow. 
That is not presently available. 
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